Anomalous drift of spiral waves in heterogeneous excitable media 
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We study the drift of spiral waves in a simple model of heterogeneous excitable medium, having 
gradients in local excitability or cellular coupling. For the first time, we report the anomalous drift 
of spiral waves towards regions having higher excitability, in contrast to all earlier observations in 
reaction-diffusion models of excitable media. Such anomalous drift can promote the onset of complex 
spatio-temporal patterns, e.g., those responsible for life-threatening arrhythmias in the heart. 
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Spatial patterns of activity, in particular spiral waves, 
are observed in a broad class of physical, biological and 
chemical excitable systems 1]. One of the most impor- 
tant contexts in which spiral waves occur is that of elec- 
trical activity in the heart, where they can act as local 
sources of high-frequency excitations. This disrupts the 
rhythmic pumping action of the heart, leading to irregu- 
larities known as arrhythmias . Understanding the dy- 
namics of spiral waves may potentially result in improved 
methods for controlling such arrhythmias ^3,, 4,, ^, J,, ^] . 
Spiral wave dynamics, primarily characterised by the mo- 
tion of its core (i.e., the trajectory of the spiral wave tip, 
defined to be a phase singularity) can be either stationary 
rotation, or, evolving with time as in the case of mean- 
dering and drift [9j- Drift, which has a significant linear 
translational component, is a possible underlyingmecha- 
nism for polymorphic ventricular tachycardia [1,1^ . This 
arrhythmia, which is characterised by an aperiodic elec- 
trocardiogram, can be a precursor of fully disordered ac- 
tivity that characterizes potentially fatal ventricular fib- 
rillation [To| . Therefore, understanding the mechanisms 
leading to spiral wave drift is not only a problem of cen- 
tral interest for physics of excitable media, but also has 
potential clinical significance. 

One of the most important causes of spiral drift is the 
heterogeneous nature of the excitable medium. This was 
first predicted in cellular automata models with step-like 
or discontinuous inhomogeneity (Tl| , which was later con- 
firmed by experiments |12l . |l3[ . Subsequently, drift has 
been shown to be induced by a smooth gradient of ex- 
citability in both simple and bi olog ically realistic ionic 
models of cardiac tissue [H, Theoretical argu- 
ments indicate that the direction of the transverse com- 
ponent of the spiral drift (i.e., orthogonal to the gradi- 
ent) depends on model parameters On the other 
hand, the longitudinal component is always directed to- 
wards the region with lower excitability (corresponding 
to longer periods of spiral rotation) ,16.] . This phe- 
nomenon has been seen in a variety of excitable media 
models of different complexity [3, 0, EBIl . However, till 
date there is no satisfactory understanding of the rea- 
sons behind the spiral wave drift towards regions with 
longer rotation period. Although earlier kinematic stud- 



ies suggested the possibility of drift towards region with 
higher excitability [l^, it has never actually been ob- 
served in a model of excitable tissue. The occurrence 
of drift in the direction of shorter period (correspond- 
ing to higher excitability) may have clinical significance, 
as it moves the spiral wave to a section where it rotates 
faster. This increases the likelihood of onset of additional 
wavebreaks away from the core, in the region where the 
medium is less excitable. Thus, it is a possible generation 
mechanism for "mother rotor" fibrillation (19, 20, 21j], 
characterised by a stationary persistent source of high- 
frequency excitations giving rise to turbulent activity in 
the heart. 

Electro-physiological heterogeneities in cardiac tissue 
may arise, in general, through spatial variation in the 
ionic currents of the excitable cells. There can also be 
gradients in the inter-cellular coupling as a result of the 
inhomogeneous distribution of the conductances of gap 
junctions connecting neighboring cells. In this paper, we 
use a simple model of cardiac tissue to investigate the role 
of both these types of heterogeneities in governing the 
direction of the spiral wave drift. We report the existence 
of a regime where the spiral wave core moves towards the 
region with higher excitability. This is a novel finding, 
never before observed in a model of excitable media and 
it may substantially increase the understanding of how 
heterogeneities affect spiral wave dynamics in the heart. 

A generic model of excitable media that describes the 
dynamics of trans-membrane potential V in cardiac tis- 
sue has the form 



dV/dt = V^DVV + ahon{V,g^ 
dg,/dt = F(V,g,). 



(1) 
(2) 



Here, lion is the total ionic current traveling through the 
channels on the cellular membrane, D accounts for the 
inter-cellular coupling and gi describes the dynamics of 
gating variables for the various ion channels. In this pa- 
per, we study the effects of heterogeneous distribution of 
the ionic currents and intercellular couplings. For this 
purpose, we introduce the parameters a and 7, which 
represent the spatial variation in ionic currents and con- 
duction properties (respectively) for an inhomogeneous 
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medium. Parameter a directly scales the value of the 
ionic current in Eq. [TJ while 7 scales the diffusion coeffi- 
cient as -D = Do+j{x) {Do = 1 for the rest of the paper). 
In this study, we have used the Barkley model [22j , where 
the several gating variables are aggregated into a single 
variable g that controls the slow recovery dynamics of the 
medium with F(V,g) = V—g. The nonlinear dependence 
of the ionic current on the fast variable V is represented 
by the cubic function /io„ = [V{1-V){V - {{g + b)/a)]/e, 
where a and b are parameters governing the local ki- 
netics and e is the relative time scale between the lo- 
cal dynamics of V and g. The spatial heterogeneity of 
local excitability and cellular coupling are assumed to 
have linear functional form, viz., a{x) = ao + Aa x and 
7(x) = 70 -l- A7 X. The variable x represents the spatial 
position along the principal direction of the inhomogene- 
ity gradient, the origin being considered to be the initial 
position of the spiral wave tip. At this point, a — ao, 
7 = 7o, and Aa, A7 measure their rate of change along 
the gradient. For all the figures in this paper, we have 
used tto = 1-15 and 70 = 1.3. 

The two dimensional system is discretized on a square 
spatial grid of size L x L (L = 200 for the figures shown 
here). The values of space step Aa; and time step At 
used are 0.5 and 0.005 respectively. A sample of simu- 
lations have been repeated for Aa; = 0.25 to verify nu- 
merical accuracy. The model equations are solved using 
forward Euler scheme with a standard 5-point stencil for 
the Laplacian. No-flux boundary conditions are imple- 
mented at the edges of the simulation domain. The initial 
condition for all simulations is a stable non-meandering 
spiral, the spiral tip being at the centre of the simulation 
domain. 

Increasing either the ionic current (via a) or inter- 
cellular coupling (via 7) results in increasing the ex- 
citability of the medium. Thus, to investigate the role 
of heterogeneity in spiral drift, we have considered spa- 
tial gradients in a or 7 individually (keeping the other 
parameter constant). After extensive numerical simula- 
tions that scan over the (a, b) parameter space of the 
Barkley model, we have found that it is indeed possible 
to observe anomalous drift of the spiral, i.e., drift towards 
regions with higher excitability (shorter periods). An ex- 
ample of such anomalous drift is shown in Fig. [1] (A,C). 
For comparison, in Fig.[T](B,D) we show the normal drift 
of the spiral towards regions of lower excitability. This 
is seen for a set of (a, b) values which is farther from the 
boundary with the sub-excitable region of the Barkley 
model [23] than the (a, 6) parameter set for which anoma- 
lous drift is observed in Fig. [1] (A,C). 

To analyse the genesis of anomalous drift, we first look 
at how the parameters 7 and a affect the spiral wave 
in an homogeneous medium. As 7 is only a scaling fac- 
tor for the diffusion coefficient, the period of the spiral 
wave does not depend on it. Further, scaling arguments 
suggest that the spiral wavelength increases as a square 
root of 7. Thus, for normal drift in the presence of a 
gradient in 7, the spiral moves towards the shorter wave- 




FIG. 1: Anomalous and normal drift of spiral wave. 

Pseudocolor images of spiral wave at the instant when the 
gradient in local excitability a (top row) and cellular cou- 
pling (bottom row) is applied. (A, C) Anomalous drift in the 
direction of higher excitability, i.e., towards increasing values 
of a (A) or 7 (C), the direction being shown by solid arrows. 
(B, D) Normal drift in a (B) or 7 (D) gradient. Parameter 
values are a = 0.82, b = 0.13 (for A, C) and a = 1.02, b = 0.15 
(for B, D). The gradients applied are (A, B) Aa = 0.0025, 
A7 = 0, and (C, D) Aa = 0, A7 = 0.020. In all cases, the 
gradient is along the vertical direction, with a or 7 increasing 
from top to bottom. In (B,D) the region around the core is 
magnified to make the wavelength of the spiral comparable to 
that in (A,C). 

length region, while for anomalous drift, it is directed 
towards longer wavelengths. Fig. [2] shows the variation 
of the spiral period and wavelength as a function of the 
parameter a, both of which decrease as a increases (2^ . 
From these results we can infer that, for normal drift in 
the presence of a gradient, the period and wavelength 
of the spiral increase as the core moves towards lower a 
regions. In contrast, we see a decrease in the period and 
wavelength in the case of anomalous drift towards regions 
having higher values of a. 

Next, we study the effect of the magnitude of spatial 
gradient in a or 7 on the velocity of spiral drift. Fig. [3] 
shows the longitudinal component of the drift velocity, 
Vl, i.e., along the gradient, as a function of the spatial 
variation in a or 7. Note that, positive Vl corresponds 
to anomalous, while, negative Vl corresponds to normal 
drift of the spiral wave. Fig. [3] shows that, for normal 
drift, increasing either of the gradients results in a mono- 
tonic increase of Vl- However, in the case of anomalous 
drift as a result of a gradient, we see a non-monotonic 
behavior in Vl, which first increases but then decreases 
and becomes negative (Fig. [31 a). Thus, the anomalous 
drift of the spiral towards higher excitability in a gradi- 
ent is seen only for small Aa. For higher Aa, there is a 
reversal of direction and the spiral exhibits normal drift. 
On the other hand. Fig. [3] (b) shows that for a gradient 
in 7, the anomalous drift is observed for the entire range 
of A7 that is investigated. 
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(a) 



■ a=0.82,b=0.13 
• a=1.02,b=0.15- 



(b) 



■ a=0.82,b=0.13 
• a=1.02,b=0.15 
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FIG. 2: The variation of spiral period and wavelength 
as a function of the parameter a. The symbols "1" and 
"2" correspond to the values of a in the region around the 
initial and final positions (respectively) of the spiral waves in 
Fig. [1] with the same sets of Barkley model parameters being 
used. The solid and broken arrows represent the directions 
of anomalous and normal drift, respectively, in presence of a 
gradient in a. 




FIG. 4: (a) Non-monotonic variation of the longitudinal com- 
ponent of drift velocity, Vl , as a function of the model param- 
eter a. The two curves correspond to media having a constant 
gradient in excitability (circles: Aa — 0.0025, A7 = 0) and 
cellular coupling (squares: Aa — 0, A7 — 0.020). For both, 
the parameter b = 0.13. (b) The contribution to Vl from the 
different components in the diffusion term as a function of the 
gradient in cellular coupling, V7. The circles and diamonds 
correspond to the linear and second-order contributions, and, 
squares correspond to the complete Laplacian term, respec- 
tively. All data points shown are for a — 0.82, b = 0.13. 

ity, but becomes normal at higher excitabihty. However, 
for the gradient in a, the anomalous drift occurs only 
over an intermediate range of a. For lower and higher 
excitabilities, the drift becomes normal. 

The mechanism of anomalous drift in a gradient re- 
mains unclear. However, we can understand the anoma- 
lous drift for the cellular coupling (7) gradient, by relat- 
ing it to other drift phenomena in excitable media. Note 
that, the Laplacian term in Eq. ([l} can be expanded as, 

yj{x)D\7V = (Do + 7(x))VV + d"f/dx dV/dx. (3) 



FIG. 3: Drift velocity depends on the gradient in pa- 
rameters a and 7. (a) Non-monotonic variation (solid 
curve) of the longitudinal component of spiral wave drift ve- 
locity Vl as a function of the gradient in local excitability, Aa, 
for a model system with parameters a = 0.82, b = 0.13. Posi- 
tive values of Vl indicate anomalous drift. For a different set 
of parameters (o = 1.02, b = 0.15), normal drift is seen for the 
entire range of gradients used (broken curve), (b) Variation 
of Vl with the gradient in cellular coupling, A7. Solid and 
broken curves represent the anomalous and normal drift seen 
for the two parameter sets mentioned earlier (respectively), 
and are observed throughout the range of gradients used. 

We have also studied the effect of the local kinetics on 
anomalous drift by varying the Barkley model parame- 
ter a (Fig. m a). Increasing a (keeping 6 fixed) decreases 
the activation threshold of the medium, and thus makes 
the system more excitable. We observe that for both a 
and 7 gradients, the variation of as a function of a 
is non-monotonic. For the cellular coupling (7) gradient, 
the presence of anomalous regime clearly correlates with 
excitability. The drift is anomalous at lower excitabil- 



Therefore, the heterogeneous cellular coupling 7(2;) con- 
tributes to both the gradient {dj/dx dV/dx), as well 
as, second order {■^{x)V'^V) terms. The relative contri- 
butions of these terms to the longitudinal component of 
drift velocity is shown in Fig. [5] (b). We observe that 
the principal effect on Vl is due to the d'y/dx dV/dx 
term, while ■y{x)V^V accounts only for about 10% of the 
observed drift. This allows us to propose the following 
explanation for anomalous drift in the presence of a gra- 
dient in 7. If we do not consider the 'y{x)V'^V term in 
the Laplacian, the spatial operators in Eq. ^ are seen 
to be identical to those in equations describing drift of a 
spiral wave in the presence of an electric field [25i] . This 
latter, in turn, is similar to the Laplacian describing the 
drift of radially symmetric filaments of a scroll ring in 
three-dimensional excitable media [1^, [13] . As shown in 
Refs. [H, , the drifts observed in these two kinds of 
systems are induced by the same instabilities. We see 
from Fig. [5] (b) that the gradient d'y/dx dV/dx, which 
determines the drift in an electric field and that of scroll 
wave filaments, also determines the drift as a result of 
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7 gradient. Therefore, we infer that the anomafous drift 
direction (towards higher excitabiHty) observed by us is 
also a result of the same long wavelength instabilities de- 
termining the drift of scroll wave filaments. This suggests 
that the occurrence of scroll expansion in 3-D implies the 
existence of anomalous drift in 7 gradient in 2-D. Con- 
versely, observation of anomalous drift might suggest pa- 
rameter regions where scroll wave expansion is possible. 

In this paper, we have explicitly demonstrated the pos- 
sibility of spiral waves to drift towards region of higher 
excitability in a simple model of heterogeneous excitable 
medium. Most of the detailed ionic models for cardiac 
tissue have the same form as Eq. ([ij and, therefore, our 
analysis can be easily extended to biologically realistic 
models, such as LRl or TNNP [sol, ^SlJ. It might be 
possible to infer the parameter range in realistic models 
where anomalous drift may occur by using the relation 
between the cellular coupling gradient induced drift and 
scroll ring expansion. Note that, the latter phenomenon 
has recently been seen in the LRl model [27[ . 

Spiral waves are not only relevant for cardiac tissue, 
but are also observed in many different excitable media. 
Thus, it may be possible to relate our observation with re- 
sults of kinematic studies llSl and models of cyclic cataly- 
sis in replicating entities (32| , which predict drift towards 
region with shorter periods. From a clinical perspec- 



tive, anomalous drift is important as it clearly promotes 
arrhythmia and may result in fibrillation by promoting 
wave-breaks away from the spiral core. Spiral drift in 
the presence of a cellular coupling gradient maybe a key 
factor giving rise to abnormal wave activity in regions 
of the heart where conductivity changes abruptly, e.g., 
at Purkinje-muscle cell junctions or in an infarct border 
zone [S^. It can also be studied experimentally and nu- 
merically in many model systems, such as, heterogeneous 
mono-layers of neonatal rat cardiomyocytes • 

To conclude, we have observed that spiral waves in 
heterogeneous excitable media can drift towards regions 
having higher excitability. Such anomalous drift occurs 
either in media having intermediate to low excitability 
when the heterogeneity is a gradient in ionic current, or, 
in media with low excitability for a gradient in the cellu- 
lar coupling. Further, it appears to be related to regimes 
where expansion of 3-dimensional scroll wave filaments is 
observed. Anomalous drift of spiral waves may increase 
the likelihood of the onset of complex spatio-temporal 
patterns in excitable medium, e.g., turbulent electrical 
activity in the heart. 
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